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The spectral properties of interacting strongly chaotic systems are investigated for growing inter¬ 
action strength. A very sensitive transition from Poisson statistics to that of random matrix theory 
is found. We introduce a new random matrix ensemble modeling this dynamical symmetry breaking 
transition which turns out to be universal and depends on a single scaling parameter only. Coupled 
kicked rotors, a dynamical systems paradigm for such transitions, are compared with this ensemble 
and excellent agreement is found for the nearest-neighbor-spacing distribution. It turns out that 
this transition is described quite accurately using perturbation theory. 


Quantization of fully chaotic systems is known to lead 
to the spectral fluctuations of random matrix theory 
(RMT) [5 and to exhibit energy level repulsion [2] . More 
generally, even non-integrable models without apparent 
classical limits, such as spin systems, can also show such 
features B SI- In contrast, integrable systems gener¬ 
ally follow Poisson statistics, which are devoid of level 
repulsion [S]. It is also well understood that combining 
spectra of different irreducible representations tends to¬ 
ward Poisson statistics in the limit of superposing many 
sequences mm- An instance where this occurs is for the 
spectra of two separable, but individually chaotic sys¬ 
tems. Whereas each subsystem possesses RMT fluctua¬ 
tions, the full spectrum tends to Poisson fluctuations in 
the large dimensionality limit [^. 

This begs the question of what happens to spectral 
fluctuations if such separable, but quantum chaotic sub¬ 
systems, interact. There are many motivations for study¬ 
ing such systems. For example, they may be of direct 
physical interest, such as conduction electrons in chaotic 
quantum dots interacting through a screened Coulomb 
potential [H]. Another motivation derives from quantum 
information theory where the development of entangle¬ 
ment is of particular importance |10j . Two quantum spin 
chains with RMT spectral fluctuations coupled in a lad¬ 
der configuration is a many-body situation where such 
transitions are possible as well. 

Typically, the interaction between two subsystems 
leads to entanglement and paves the way for RMT fluc¬ 
tuations of the combined system. This Poisson-to-RMT 
transition can be viewed as a dynamical symmetry break¬ 
ing in analogy to fundamental symmetry breaking; the 
first exact RMT solution to a symmetry breaking prob¬ 
lem involved time-reversal invariance m- For modeling 
a particular dynamical system, it is important to connect 
dynamical system parameters with the abstract transi- 
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tion parameter built into an RMT transition ensemble, 
and to confirm that this correctly describes the dynam¬ 
ical system’s statistical properties. For example, partial 
transport barriers often arise classically and an RMT 
ensemble can correctly model their effects, but only if 
they have been accounted for their respective fluxes, rel¬ 
ative phase space volumes, and tendency to couple lo¬ 
cally [nma. 

There are a couple of other possible cases of a Poisson- 
to-RMT transition, a metal-insulator transition where 
states transition from localized to extended and 

the perturbation of an integrable dynamical system m 
II5j that renders it chaotic for sufficient perturbation 
strength. In neither of these possibilities is there a sim¬ 
ple globally coupled Poisson-to-RMT ensemble. Various 
complications, such as the metal-insulator transition, the 
KAM theorem regarding the survivability of tori, and 
partial transport barriers, all prevent this. 

In this paper we address how strong the interactions 
must be between the subsystems in order to recover the 
spectral fluctuations of fully chaotic systems, and how 
the Poisson-to-RMT transition occurs as the interaction 
magnitude varies from non-interacting to strongly inter¬ 
acting particles cases. We show that the interplay be¬ 
tween the coupling and an effective Planck constant gives 
rise to a dimensionless scaling transition parameter. This 
is obtained in much the same way as found in the con¬ 
text of a global symmetry breaking in RMT nmM]. 
There, arbitrarily small couplings, but uniform every¬ 
where, lead discontinuously to level repulsion [21j . If 
scaled properly, the entire transition is universal and pre¬ 
dicted well by a perturbation theory. 

Interacting systems .— To study the Poisson-to-RMT 
transition, consider the unitary Floquet operator, or 
kicked version, of generic bipartite systems described by 

U={Ux® C/ 2 ) C /12 (1) 

where the Uj are the subsystem Floquet operators and 
C /12 describes the interaction. An example is the Hamil¬ 
tonian 

H = +P 2 ) + [U('Zi) + k2((72) + bVi2{qi,q2)\ St (2) 



2 


where 5t = J2’^=-oo 5{t — n) is a periodic train of kicks, 
and a unit time has been chosen as the kicking period. 
The propagator connecting states separated by one kick 
follows from Uj = exp[—ip'j/{2h)]exp{—iVj/h) (j = 1,2) 
and Ui 2 = exp{—ibVi 2 /h). The classcial limit is a 4- 
dimensional symplectic map {qj,pj) i— {q'j,p'j) = {qj + 
p'j,Pj — dVjjdqj — b dVi 2 /dqj), connecting the state of 
the system immediately prior to consecutive kicks. 

Perhaps the simplest paradigm for our purposes is that 
of coupled kicked rotors [22]. They have been realized 
in experiments on cold atoms |53|, can be made strongly 
chaotic, and the interaction strength continuously varied. 
The most elementary case is for two interacting rotors 
with the single particle potentials Vj = Kj cos(2Trqj)/ 
and interaction 

Vi 2 = ^ cos[27r(gi + 92 )]- (3) 

The unit periodicity in the angle variables qj is extended 
here to the momenta pj so that the phase space is a 4- 
dimensional torus. If the kicking strengths, {Ki,K 2 }, 
here {9,10} respectively, are each chosen sufficiently 
large, the individual maps are strongly chaotic with a 
Lyapunov exponent of « ln(iT,/2) [24|. The interaction 
strength is tuned by the parameter b. 

The quantum mechanics on a torus phase space of 
a single rotor gives rise to a finite Hilbert space of di¬ 
mension N. The effective Planck constant is h = 1/N. 
Thus Ui and U 2 are N dimensional unitary operators 
on their respective spaces Hi and H 2 whereas the in¬ 
teraction U 12 is a iV^ dimensional unitary operator on 
the tensor product space H^ 0 H 2 ■ The quantized 4- 
dimensional map is given by Eq. Q. Such coupled quan¬ 
tum maps have been studied in different contexts [53 [2^ 
where more details can be found. Here we are inter¬ 
ested in the statistics of the eigenphases Pn defined by 
^IV’n) = The boundary conditions are chosen 

to break both parity and time-reversal symmetries. Note 
that coupled kicked rotors with a different interaction 
term have also been studied on the cylinder m- 

RMT transition ensemble. — The first task is to con¬ 
struct the RMT transition ensemble associated with 
Eq. 0, 

Urmt (e) = ® ^^cue) Ui2ie) (4) 

where the tensor product is taken of two independently 
chosen circular unitary ensembles (CUE) of dimension N 
representing the individual strongly chaotic rotors, and 
Ci 2 (e) a diagonal unitary matrix in the resulting N'^ di¬ 
mensional space representing the coupling. Its diagonal 
elements are taken as exp(27rie^„^„2) where ^nin 2 are in¬ 
dependent random variables that are uniform on the in¬ 
terval (—1/2,1/2], e is a real number, and 1 < ni, 712 < 
label subsystem bases. Any special structure of the ro¬ 
tor’s C /12 is ignored because it turns out to be irrelevant 
for this study. A limiting case of this ensemble has been 
studied previously uni, wherein the entangling power of 


Z^RMT(e = 1) was found analytically. However, its spec¬ 
tral features for general e are not yet explored. 

For simplicity, the only statistical quantity considered 
in this paper is the nearest neighbor spacing (NNS) dis¬ 
tribution P{s) of the eigenphases (/?„. For e = 0, C/rmt 
is the direct (tensor) product of two independent CUE 
matrices. In this case P{s) has been shown analytically 
as A^ —>■ 00 to approach the Poissonian result of exp(—s), 
where s is the spacing in terms of unit mean [H|. The 
case e = 1 represents strong coupling and CUE behavior. 

Transition of the NNS distribution .— For non¬ 
interacting rotors, the Poisson result is also expected to 
hold because the separability leads to spectral sequences 
that are composed of superpositions of sequences. For 
interacting rotors, the relation between e of Eq. and 
the b of the coupled rotors is needed. As explained below. 



for small values of e and Nb, and A^ ^ 1. In practice 
this approximation is good for the entire transition even 
if N is only moderately large. 

Figure[^shows a comparison of the results for the RMT 
transition ensemble and two coupled rotors using Eq. (HI- 
This demonstrates that the ensemble captures the tran¬ 
sition in statistics from Poisson-like towards global CUE 
perfectly well. CUE statistics already arise for small val¬ 
ues of &, which reflects the sensitivity of the transition. 

Universal scaling .— The second task is to identify a 
universal scaling parameter, if it exists, and test whether 
it captures the proper transition scale. In other instances 
of symmetry breaking miiiiiiHi US], it emerged using 



FIG. 1. (Color online) Nearest neighbor level-spacing distri¬ 
bution P{s) for two coupled rotors (circles), and the RMT 
model (triangles), (i) For & = 0, it is Poissonian (red line) 
as is the e = 0 RMT ensemble, (ii) For b = 0.0019, the 
spacing distribution is intermediate and agrees well with the 
corresponding RMT ensemble with e = 0.012 and the result 
of a perturbation theory (blue line) for A = 0.1158. (in) For 
b — 0.008, the spacing distribution agrees well with both the 
RMT ensemble (e = 0.05) and the CUE. For all cases the 
dimensionality of the unitary matrix is N'^ = 10^^. 
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perturbation theory where it was given by A = jD'^, 

with the mean square off-diagonal matrix element in 
the unperturbed system’s diagonal representation and D 
its mean level spacing. If A ^ r, where r is the effective 
range of the statistics (for the NNS, r ~ 1 ,) the transi¬ 
tion is nearing completion. It must be emphasized that 
the interaction is assumed to be generic, has no special 
symmetries and is entangling. Ideally is calculated 
only from those matrix elements that are responsible for 
the spectral transition, but averaging over a global set 
of all matrix elements gives the identical result if there 
is no secular structure, e.g. like bandedness or special 
correlations [ 2 H] . 

For an dimensional unitary matrix, necessarily 
D = 2'k/N'^. That leaves v'^ to calculate. The un¬ 
coupled part of the transition ensemble is diagonalized 
by the direct product of two independent N dimen¬ 
sional unitary matrices, say u and w. An off-diagonal 
element of the interaction operator in the unperturbed 
basis is Zki-k'u = m,n'^kmU*k,^winWp^{U 12 )mn, where 
{k,l) 7 ^ {k',l'), and {Ui 2 )mn = {mn\Ui 2 \mn) is a diag¬ 
onal matrix element due to the interaction. All indices 
refer to subsystems and therefore lie in [1,A^]. The off- 
diagonal elements are restricted to k ^ k' and I ^ I', as 
the energies of neighboring states after the tensor prod¬ 
uct must be different in both. 

The required average is that of \zki-k'i'\‘^, as u and w 
are independently chosen according to the Haar mea¬ 
sure on U(N) , the group of unitary matrices in the 
subsystem spaces. This is done exactly using known 
results for the average of which is 

{Smm' — ^/N)/{N'^ — l) whenfc 7 ^ k' (for example see P5]l. 
and leads to (for any N) 


- 47r2(Ar2 _ 1)2^ 


/ 

1 

^ 1 

r 




||t/Wf + ||[/(2)|| 



( 6 ) 


Here U^P = J2kiUi2)ik, ujp = J2iiUi2)tk are partially 
traced (still diagonal) interaction operators, which are in 
general not unitary, and |i ATp = Tr(XvAl) is the Hilbert- 
Schmidt norm. Performing the exact ensemble average 
over t/i 2 (e) gives 


A = 


N* 



sin^(7re) 


iV2g2 

12 


(7) 


where for the last approximation A^ ^ 1 and e <C 1 
have been used. Note that related averages involving 
Haar measures over product spaces have been obtained 
recently [50] . 

It is important to recognize that Eq. in the large- 
N limit, applies to individual chaotic dynamical systems, 
and hence the coupled rotors. The eigenbasis of uncou¬ 
pled rotors behaves like that of a CUE member, and in 
the large-A^ limit has the same statistical behavior as the 


ensemble. Thus, applying Eq. ([^ to the coupled rotors 
gives 


A=2^{l-jPNb/2.)) 


327r4 ’ 


( 8 ) 


where Nb ^ 1 in the right form. Using the approxima¬ 
tions in Eqs. (^) generates Eq. ([^. It can be checked 
that this A correctly vanishes if V 12 is a function of either 
coordinate alone, and thereby not entangling. However, 
it does not vanish for general separable potentials, and 
hence the assumption that the interaction is entangling 
is necessary. 

There are two significant consequences of Eqs. ffl- 
The first is that there is a universal transition governed 
by the scaling parameter A oc oc b^N'^. Indeed, this 
is well verified in Fig. Varying b and N, but holding 
b'^N'^ fixed generates precisely the same statistical fluctu¬ 
ations. Given that RMT ensemble statistics match those 
of the coupled rotors. Fig. the transition of Eq. Q is 
the same for fixed this was also verified, but left 

out of Fig. [^ to avoid clutter. Secondly, the transition 
is very sensitive to the interaction strength. The further 
a system is in the short wavelength limit, i.e. A —> 0 , 
the smaller the interaction strength needed to drive the 
transition to completion. This is seen in the small val¬ 
ues of b,e chosen in Figs. |1|2| In fact, the transition is 
discontinuous in the interaction strength as A^ —^ 00 m- 

For Hamiltonians of the Eq. ([^-type and perturbative 
interactions (second order in nV), Eq. simplifies to 

A = [{VP)^2 + (Ui2)?2 - ((^12)?)2 - ((Vl2)^)l] . 

(9) 

Here {X)j = Tr^V/V, {X)i 2 = Tt:i 2 X/N'^. For the cou¬ 
pled rotors, with the interaction as in Eq. ([^, the last 
three terms vanish in the large N limit as the integral over 



FIG. 2. (Color online) Universal scaling of the level-spacing 
distribution P{s) for two coupled rotors. Three results (tri¬ 
angles, squares, circles) are shown for {b,N) = (0.0025,100), 
(0.01,50), and (0.04,25). They all correspond to the same 
transition parameter A « 0.2. The solid curve comes from 
a perturbation theory described in the text and the dashed 
curves are the Poisson and CUE distribution, respectively. 
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qi or (72 is zero. Under such conditions, A « {^^ 12 )^ 

where (X) = Jq Jq dqidq 2 X(qi, q 2 ) and X is a classi¬ 
cal function of the two coordinates. The interaction is 
assumed to be diagonal in the position basis for the con¬ 
tinuum approximation. 

Perturbation theory .— A perturbative theory can be 
derived based on previously developed techniques na 
[3l]. The first step is to consider a very large dimensional 
case and the perturbation expansion for the eigenvalues 
to second order. The first order term is diagonal and 
shifts the energies around randomly, and thus does not al¬ 
ter the Poisson nature of the spectral statistics. Amongst 
the second order terms, terms connecting lower energies 
push a level up and those above push down. Overall, a 
level’s motion fluctuates, but remains in the same neigh¬ 
borhood. However, when calculating the NNS spacing, 
there is a common term, which enters opposite in sign 
for the two levels, and always pushes them apart. To a 
rough approximation, the bulk of the perturbation terms 
keep the mean spacing constant, but the exceptional term 
introduces level repulsion. 

Thus if AE° is the spacing between two neighbor¬ 
ing levels at (i,i -I- 1) before perturbation, it becomes 
AEi « AE^ + 2|Uy+ip/AU°, ignoring the effect of 
other levels. Here is an interaction matrix element, 

whose real and imaginary parts are complex Gaussian 
random variables with a scale determined by the vari¬ 
ance introduced above. This is appropriate given the 
normal fluctuations of the real and imaginary parts of the 
GUE matrix elements. Scaling the spacings and matrix 
elements, r = AEi/D,so = AE^/D and = v'^w, 

introduces A as above, r « sq + ‘2Aw/so, with w and sq 
distributed exponentially, the latter due to the Poisso- 
nian nature of the unperturbed spectrum. Close lying or 
nearly degenerate levels are frequent for an unperturbed 
Poissonian spectrum. The resulting divergences are regu- 



FIG. 3. (Color online) Variance of the nearest neighbor spac¬ 
ing as function of the transition parameter. Shown are results 
for the coupled rotors (circles) and the RMT model (trian¬ 
gles), in both cases for N = 50. The perturbation theory 
result is shown as solid line, Eq. (131, while the dashed line 
refers to a 2 x 2 transition ensemble result 1321. 


larized using degenerate perturbation theory, which gives 

p{r) = J dso dw exp(—So — w)S ~ ■\/so + dAic^ 

( 10 ) 

The result of integration is 


p{r) 


Va 




/r-2A 

V 2Va 


^2 

-I- e“ JA £)_|_ 



( 11 ) 


where D+{x) is the Dawson function dt . The 

effect of all the other levels is accounted for by recom¬ 
pressing the spectrum to unit mean spacing using 


(r) 


1-b A 



- e-^Ei (A) 


( 12 ) 


where Ei(a;) is the exponential integral. Thus, the 
rescaled NNS distribution is P(s) = (r)p(s(r)). Also, 
the variance of the NNS is 


a 


2 


2 + 4A 


(13) 


and is a useful parameter to monitor the whole transition. 
Figure [^compares the NNS obtained using Eqs. ( [II| ) and 
© with those of the coupled rotors, and is an excellent 
approximation. Figure compares the variance across 
the whole transition. While the perturbation theory re¬ 
sult is slightly better than a 2 x 2 model results |32] for 
A < 1.0, the latter by design does better for larger A. 
A more complete theory for the NNS covering the whole 
transition is left for future work. 

Summary and outlook .— In summary, this paper ad¬ 
dresses an important, rather general question of what 
happens when two chaotic systems, or systems with ran¬ 
dom matrix fluctuations interact. The resulting entan¬ 
glement leads to a universal transition and rapid recovery 
of global RMT fluctuations. The transition is very well 
captured by a natural random matrix ensemble, Eq. (|^ . 
A universal scaling parameter is derived in terms of the 
interaction strength and separately the abstract coupling 
parameter of the RMT ensemble, thus establishing a 
quantitative relationship between the two objects. Fi¬ 
nally, perturbation theory with a simplified global Pois¬ 
son plus RMT ensemble reproduces the NNS transition 
in fluctuations very well. In spite of its application to a 
particular dynamical system, the coupled kicked rotors, 
the methods and results are of a general nature, and can 
be expected to hold wherever there is a global coupling 
of two subsystems each exhibiting RMT fluctuations. 

Various extensions of this work are of significant in¬ 
terest. One rather immediate generalization is to sub¬ 
systems with different dimensions, say a bipartite sys¬ 
tem with differing numbers of spins in the subsystems. 
In many of the relations, the substitution —>■ N 1 N 2 

suffices to capture the new dependences. A much more 
involved change is to a greater number of subsystems. If 
the number of subsystems is L, each with a dimension¬ 
ality N, a new expression for the transition parameter 
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equation in Eq. is possible. There, the new RMT 
model analogous to Eq. Q would have L CUE matri¬ 
ces in a tensor product and a global diagonal coupling of 
the same kind. That leads to a generalization of Eq. 0: 
A = N^^{l-smc^{Tre))/[4TT^{N + l)^] « jVl. Thus 
for large N, it is expected that the recovery of RMT fluc¬ 
tuations is much faster with increasing h. However, this 


analysis assumes essentially that Ul{^) is an L-body op¬ 
erator. If it is restricted to a lower rank form, say a 
2-body operator, more sophisticated models have to be 
separately considered. Preliminary results have also been 
obtained regarding the statistics of eigenfunctions, and 
entanglement in the eigenstates which develop solely due 
to the interactions. They are beyond the scope of this 
paper and left for future publication. 
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